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Abstract 


Compression  of  a  solid  material  under  plane  shock  wave  propagation/  loading 
can  provide  information  related  to  change  in  density  and  shear  strength  with 
pressure/  stress  of  the  material  under  inertial  confinement.  The  present  work 
surveys  and  analyzes  the  existing  shock  compression  and  high-pressure  data  on 
silicon  carbide  to  provide  pertinent  information  to  material  modelers  and  others 
for  use  in  their  work  related  to  this  material.  It  also  points  out  a  set  of 
investigations  that  needs  to  be  carried  out  to  complement  the  existing  data. 
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1.  Introduction 


Silicon  carbide  (SiC)  is  a  material  of  considerable  interest  because  of  its 
technological  applications  as  superabrasive,  high-temperature  electronic 
material,  as  well  as  blue  light  emitting  diode,  and  it  is  a  material  component  of 
an  armor  system  [1-3].  AU  SiC  structures  are  made  up  of  a  pair  of  hexagonal 
layers,  one  made  up  of  carbon  atoms  and  the  other  of  silicon  atoms.  The  two 
simplest  systems  for  SiC  are  cubic  zincblende  (P-SiC)  and  wurtzite  structure  (a- 
SiC).  Both  have  densities  of  3.21  Mg/m^.  p-SiC  transforms  to  a-SiC  above 
1873  K.  SiC  is  also  reported  to  exist  in  over  100  polytypes  tmder  the  ambient 
condition  [4].  The  different  polyt3^es  are  generated  by  the  stacking  sequence  of 
a  nxunber  (N)  of  SiC  pair  layers  in  the  cubic  [111]  direction  within  a  unit  ceU  of 
tire  cubic  (Q,  hexagonal  (H),  or  rhombohedral  (R)  structure.  Numbers  are 
attached  to  each  polytype  according  to  the  number  of  pair  layers  witiiin  the 
stack.  These  pol5^types  differ  witii  respect  to  temperature  at  generation  [5].  2H- 
SiC  and  3C-SiC  are  formed  at  1773  and  1873  K,  respectively.  4H-SiC,  6H-SiC,  and 
others  are  formed  at  and  above  2273  K.  6H  is  the  most  stable  of  these  polytypes 
at  temperatures  between  2473  and  2873  K.  Polyt5q)es  2H,  3C,  4H,  6H,  and  15R 
occur  most  frequentiy.  As  yet,  a  complete  set  of  values  of  all  the  elastic  constants 
of  either  form  of  SiC  is  not  available. 

Several  investigations  deal  with  the  hydrostatic  compression  [6-9]  of  SiCs.  In 
addition,  over  the  past  four  decades  a  large  number  of  shock  wave  compression 
experiments  have  been  conducted  on  polycrystalline  SiCs  to  determine  their 
compressibility  and  shear  strength  retained  under  plane  shock  wave 
compression  [10-20].  This  list  excludes  shock  wave  investigation  studies  dealing 
with  spallation  in  SiC  because  Dandekar  and  Barti<owski  [21]  recently  covered 
this  subject.  The  t5q>e  of  SiC  used  in  earlier  shock  compression  investigations 
was  not  always  unambiguously  identified  except  in  terms  of  densities  and  in 
some  cases  chemical  compositions.  But,  if  a  description  of  processes  employed  to 
produce  SiC  was  provided  and  then  based  on  the  stability  of  the  various 
polytypes,  one  may  infer  structure  and  polytype  from  such  a  description,  in 
recent  years,  the  t5q>e  of  SiC  used  in  experiments  is  generally  identified,  expUcitiy 
making  comparisons  of  the  results  obtained  by  the  various  investigators  easier 
and  more  meaningful.  Such  explicit  identification  also  provides  a  margin  of 
error  in  the  values  of  the  various  parameters  Kke  compressibility  of  the  material 
and  shear  strength  retained  by  a  specific  type  of  SiC  xmder  plane  shock  wave 
propagation.  This  work  summarizes  the  resiilts  of  the  existing  hydrostatic 
compression  data  and  shock  wave  compression  data  on  silicon  carbide  in  the 
open  literature. 
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2.  Compression  of  SiC 


2.1  Hydrostatic  Compression 

The  results  of  hydrostatic  compressions  of  3C-SiC,  6H-SiC,  and  15R-SiC, 
determined  by  high-pressure  x-ray  diffraction  experiments,  are  reported  by 
Strossner  et  al.  [6],  Aleksandrov  et  al.  [7],  Bassett  et  al.  [8],  and  Yoshida  et  al.  [9]. 
These  investigators  used  diamond-anvil  cells  to  generate  high  pressure  in  SiCs. 
The  pressure  attained  in  the  material  was  estimated  from  the  calibration  data  of 
Mao  et  al.  [22],  or  from  volumetric  compression  of  gold  using  the  equation  of 
state  of  Heinz  and  Jeanloz  [23].  Whereas  Aleksandrov  et  al.  [7]  confined  SiC  in  a 
compressed  helium  medium  that  produced  a  hydrostatic  pressure  environment 
estimated  to  be  over  30  GPa,  Bassett  et  al.  [8],  Strossner  et  al.  [6]  and  Yoshida  et 
al.  [9]  confined  their  materials  in  a  mixture  of  sodium  chloride  and  gold,  and 
methanol,  ethanol,  and  water,  respectively,  to  produce  a  hydrostatic  pressure 
environment  in  their  samples.  The  results  of  these  investigations  on  crystaUine 
SiC  are  summarized  in  Table  1. 

Table  1.  Bulk  modulus  (Bo),  pressure  derivative  of  the  bulk  modulus  (B'o),  and  maximum 
pressure  (Pmax)  to  which  experiments  were  carried  out  m  SiC  in  various  static 
pressure  investigations. 


Type 

Bo 

(GPa) 

B’o 

P  max 

(GPa) 

Reference 

3C 

248  ±9 

25 

[6] 

3C 

4.1  ±  0.1 

42 

m 

15R 

224  ±3 

4.3  ±0.1 

45 

.  [7]  .  _  _ 

6H 

68 

[8] 

3C 

260  ±9 

2.9  ±0.3 

105 

[9] 

6H 

260  ±9 

2.9  ±  0.3 

95 

[9] 

It  is  of  interest  to  deal  with  the  status  of  elastic  constant  measurements  in  various 
SiC's.  Lambrecht  et  al.  [24]  gave  a  complete  and  critical  review  of  experimental 
measurements  of  elastic  constants  of  SiC.  There  has  been  no  direct 
determination  of  the  full  set  of  elastic  constants  for  either  3C-SiC  or  hexagonal 
polytypes  of  SiC,  as  single  crystals  of  3C-SiC  and  6H-SiC  have  been  generally 
unavailable.  Thus,  estimates  of  elastic  constants  of  a  polycrystaUine  SiC  must  be 
made  on  some  assumption.  Hofmann  et  al.  [25]  used  the  bond-charge  model  in 
conjunction  with  the  measured  Raman  spectra,  infrared  absorption,  and 
reflectivity  measurements  and  luminescence  studies  on  3C,  2H,  4H,  6H,  lOH, 
15R,  and  21R  SiC  to  obtain  the  elastic  constants  of  various  polytypes  of  SiC.  They 
calculated  the  values  of  elastic  constants  obtained  from  the  limited  sound  wave 
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velocity  measurements,  resonance  experiments,  and  optical  and  phonon 
measurements.  Their  results  indicate  that  the  values  of  elastic  constants  of  3C, 
2H,  4H,  and  6H-SiC  are  fairly  close  to  one  another.  The  estimated  values  of  bxilk 
modulus  of  poly  crystalline  3C-SiC  range  between  225  and  247  GPa  [24].  The 
estimated  value  of  the  bulk  modulus  for  polycrystalline  6H-SiC  based  on 
experimental  values  of  single  crystal  elastic  constants  Cn,  C33,  and  C12  measured 
by  Arlt  and  Schodder  [26],  and  the  estimated  value  of  C13  by  Hofmann  et  al.  [25] 
is  273  GPa.  It  is  difficult  to  calculate  margin  of  errors  in  these  calculated  values 
of  bulk  moduli  for  these  two  pol)d3^es  of  SiC.  The  measured  value  of  the  bulk 
modulus  of  polycrystaUine  SiC-B,  which  is  primarily  6H-SiC,  from  the  ultrasonic 
measurements  reported  by  Dandekar  and  Bartkowski  [21]  is  221  ±  2  GPa  (Table 
2).  This  table  also  gives  die  properties  of  two  other  types  of  polycrystalline  SiC. 
Holmquist  et  al.  [27]  provide  these  properties  for  other  SiCs.  The  accuracies  of 
sound  velocities  and  scatter  in  the  sound  wave  velocities  between  samples  of  KT- 
SiC  were  each  reported  to  be  within  2%  [12].  The  properties  of  SiC-B  and  SiC-N 
given  in  here  were  reported  by  Dandekar  and  Bartkowski  [21]. 


Table  2.  Properties  of  various  polycrystaUine  SiCs. 


Property 

KTSiC[12] 

SiC-N  [21] 

Grain  Size  (jim) 

— 

2-10 

1-8 

3.09 

3.215  ±  0.002 

3.227  ±0.001 

Pore  Volume  Fraction 

0.04 

0.002 

— 

Elastic  Wave  Velocity  (km/  s) 

— 

— 

— 

Longitudinal 

11.4 

12.198  ±  0.026 

12.262  ±0.001 

Shear 

7.27 

7.747  ±0.018 

7.77  4±  0.005 

Bulk 

7.71 

8.29  ±0.03 

8.354  ±  0.006 

Elastic  Modulus  (GPa) 

— 

— 

Yoimg's 

378 

448.4  ±2.1 

454.0  ±  0.6 

Shear 

163 

193.0  ±  0.9 

195.0  ±0.2 

Bulk 

184 

221.1  ±1.8 

225.2  ±  0.3 

Poisson's  Ratio 

0.157 

0.162  ±0.003 

0.164  ±  0.001 

The  values  of  Bo  for  both  3C  and  6H  SiC  differ  significantly  from  each  other  in 
these  hydrostatic  compression  investigations  (Table  1).  The  reasons  are  difficult 
to  single  out  with  any  confidence.  Most  likely,  the  differing  pressure  medium 
used  in  these  investigations  leads  to  a  departure  from  differing  amoxmts  of 
hydrostaticity  in  tibese  investigations  at  differing  pressure  levels.  The  calculation 
of  pressure  at  a  given  compression  can  easily  lead  to  the  observed  difference 
especially  in  hard,  brittle,  and  incompressible  material  hke  SiC.  For  example, 
Aleksandrov  et  al.  [7]  found  that  the  magnitude  of  pressure  at  a  given 
compression  for  3C-SiC  calculated  from  their  own  calibration  and  by  adopting 
the  one  suggested  by  Mao  et  al.  [22]  differ  by  as  much  as  4  GPa  at 
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14%  compression.  The  values  of  pressure  at  this  compression  are  40  and  44  GPa, 
resulting  in  a  10%  variance  with  calibration  method.  This  also  implies  that  the 
estimates  of  B'o  from  these  types  of  experiments  may  be  subject  to  an  even  larger 
variability.  Further,  because  Segletes  [28]  showed  that 

Tsuter  =  -1/6  +  0.5  B’o,  (1) 

the  variation  in  B'o  is  on  scale  of  F  variations,  which  are  quite  large  and  can  be 
the  source  of  even  larger  variability. 

The  compression  curve  obtained  for  SiC  from  these  studies  may  be  presented 
either  in  the  form  of  the  Birch-Mumagham  finite  strain  equation  of  state,  as  is 
customary  and  favored  by  the  investigators  in  die  field  of  static  high  pressure,  or 
the  equation  of  state  based  on  the  linear  shock  velocity  (Us)  and  particle  velocity 
(up)  relationship  observed  in  solids  in  the  field  of  shock  wave  compression. 

Birch  [29]  showed  that  the  isentropic  compression  using  the  Eulerian  definition 
of  finite  strain,  i.e.. 


is  given  by 


/=[(Vo/V)2/3-l]/2, 


(2) 


P  =  3  Bof(2 /+l)2-5  [1+ai  /+  a2p+  ....],  (3) 

where  V  is  reciprocal  of  density,  a's  are  constants,  and  P  is  the  pressure.  The 
subscript  0  indicates  the  initial/  ambient  specific  volume. 

The  equation  of  state  based  on  the  Linear  relationship  between  shock  velocity  (Us) 
and  particle  velocity  (up),  given  by 


Us  =  Co  +  s  Up 


(4) 


yields 


Ph  =  Po(Co)2ti/(1-sti)2,  (5) 

where  t]  =  1  -  (V /Vo)  and  Ph  is  the  Hugoniot  pressure. 

Under  the  Birch-Mumaghan  assumption,  a2  =  0,  Jeanloz  and  Grover  [30]  showed 
that 


ai  =  1.5(B'o  -  4)  =  (6s  -  7.5)  (6) 

and 

a2  =  215/  6  -  60s  -  3  yos  +27  s^.  (7) 


In  equation  (7),  yo  is  the  Gruneisen  parameter  at  the  ambient  condition.  Further 
Jeanloz  and  Grover  [30]  showed  that  equations  (4)  and  (6)  are  identical  for  a 
material  if  the  values  of  yo  and  B'o  lie  between  the  range  of  0.8-2.6  and  4-6, 
respectively.  It  should  be  noted  that  if  terms  involving  a2  are  neglected  and  the 
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value  of  B'o  determined  from  the  static  high-pressure  experiment  is  4,  then 
equations  (4)  and  (6)  will  give  identical  results  if  and  only  if  ttie  value  of  s  from 
shock  wave  experiments  is  determined  to  be  1.25.  In  case  this  is  not  true,  the 
compression  curves  obtained  from  the  static  and  shock  measurements  will  differ 
from  one  another.  The  significance  of  this  difference  at  a  given  compression  of  a 
material  is  application  dependent.  For  example,  the  magnitude  of  pressures  for 
3C  and  6H  SiC  at  23%  compression  obtained  from  these  two  equations  for  the 
experimental  results  reported  by  Yoshida  et  al.  [9]  are  98  and  101  GPa, 
respectively. 

Figure  1  shows  die  mean  volumetric  compression  of  6H-SiC  obtained  by  Bassett 
et  al.  [8]  and  Yoshida  et  al.  [9].  The  experimented  data  coordinates  are  not  shown 
in  tiiis  figure  to  preserve  the  clarity.  This  figure  also  shows  the  lower  and  the 
upper  hmits  of  die  volumetric  compression  of  6H-SiC  based  on  the  results  of 
Yoshida  et  al.  [9].  The  lower  and  the  upper  limits  of  the  compression  are 
obtained  using  the  values  of  Bo  and  B'o  as  251  GPa  and  2.6,  and  269  GPa  and  3.2, 
respectively  (Table  1).  This  figure  shows  that  the  compressions  of  6H-SiC 
determined  in  these  two  investigations  are  very  close  to  one  another  even  though 
the  values  of  bulk  modulus  and  its  pressure  derivatives  determined  in  tiiese  two 
investigations  are  beyond  the  range  of  experimental  uncertainties  reported  by 
these  investigators.  Bassett  et  al.  [8]  assumed  the  value  of  the  pressure  derivative 
of  the  bulk  modulus  to  be  4. 


Figure  1.  Compression  of  6H-SiC  from  static  high-pressure  x-ray  diffraction  data. 
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Finally,  Yoshida  et  al.  [9]  found  that  3C-SiC  undergoes  a  phase  transition  to  a 
rock  salt  type  structure  at  and  above  100  GPa,  with  an  accompanying  volume 
reduction  of  20.3%.  This  phase  transition  is  reversible  and  goes  back  to  the 
original  3C  phase  when  the  pressure  is  lowered  to  35  GPa.  Thus,  the  persistence 
of  the  rock-salt  phase  down  to  35  GPa  from  100  GPa  amounts  to  a  large 
hysteresis  accompanying  the  phase  transition.  They  did  not  observe  any 
indication  of  phase  transition  in  6H-SiC  to  95  GPa. 

Liu  and  Vohra  [31]  reported  the  results  of  Raman  study  in  6H-SiC  to  95  GPa  in  a 
diamond  cell  anvil.  Based  on  the  Raman  spectra  of  6H-SiC,  they  concluded  that 
the  6H-SiC  structure  is  stable  to  95  GPa.  Their  sample  became  optically 
transparent  at  this  pressure.  But  the  relative  softening  of  the  Gruneisen 
parameter  for  the  longitudinal  optical  (LO)  mode  compared  to  the  transverse 
optical  (TO)  mode  at  elevated  pressure  suggested  to  them  that  it  might  be  a 
precursor  to  a  phase  transition  in  the  range  of  100-200  GPa.  Further,  their 
observation  of  initial  increase  in  transverse  effective  charge  to  40  GPa  and  an 
observed  decline  upon  further  increase  in  pressure  suggested  that  6H-SiC  is 
more  ionic  when  compressed  to  40  GPa,  but  becomes  more  covalent  upon  further 
increase  in  the  pressure.  Whether  this  leads  6H-SiC  to  become  more 
incompressible  and/or  retain  lower  or  higher  shear  stresses  at  shock-induced 
stress  above  40  GPa  remains  to  be  ascertained  through  future  investigations. 

2.2  Hydrodynamic  Compression 

Shock  compression  of  a  material  is  affected  by  the  magnitude  of  porosity  and 
impurities  present  in  it,  which  is  not  the  case  for  in-situ  high-pressure  x-ray 
diffraction  measurements  discussed  previously.  McQueen  et  al.  [10]  and  Marsh 
[11]  published  shock  compression  data  on  SiC,  with  density  varying  between 
2.33  and  3.12  Mg/  m^.  Gust  et  al.  [12]  measured  shock  compression  of  SiC  with  a 
density  of  3.09  Mg/ m^.  The  material  used  by  Gust  et  al.  [12]  was  hot-pressed 
and  manufactured  and  marketed  by  Carborundum  Company  under  the  name 
Type  KT  SiC.  They  reported  the  impurities  in  KT  SiC.  The  material  was  a-SiC, 
but  no  identification  about  the  polytype  was  provided.  Based  on  the  similarity  of 
the  highest  density  of  SiC  used  by  McQueen  et  al.  [10],  it  is  presumed  that  their 
material  was  also  obtained  from  the  Carborundum  Company.  Sekine  and 
Kobayashi  [13, 14]  recently  measured  shock  compressions  of  6H-SiC  and  3C-SiC 
to  160  and  150  GPa,  respectively.  6H-SiC  material  used  was  manufactured  by 
Cercom  Inc.  and  is  marketed  as  SiC-B.  Sekine  and  Kobayashi  did  not  mention 
the  source  of  3C-SiC,  and  only  two  experiments  were  performed.  The  results  of 
two  experiments  on  3C-SiC  were  reported  to  be  consistent  with  those  of 
Yoshida  et  al.  [9].  The  Hugoniot  Elastic  Limit  (HEL)  of  3C-SiC  was  determined 
to  be  12  GPa.  Grady  [15]  attempted  to  measure  the  hydrodynamic 
compressibility  of  SiC  through  shock  compression  of  mixture  samples  of 
a-silicon  carbide  powder  and  copper.  The  shock  compression  data  reported  by 


6 


Sekine  and  Kobayashi  [13, 14]  complements  the  results  of  low  magnitude  shock 
compression  experiments  designed  to  measure  stress  wave  profiles  and  obtain 
the  associated  properties  of  SiC-B  and  other  a-SiC,  such  as  the  HEL,  and  the 
nature  of  inelastic  deformation,  shear  strength  retained  under  plane  shock  wave 
compression,  and  release  behavior  reported  in  references  [16-20]. 

Hydrodynamic  compressions  of  these  materials  represented  by  the  bulk 
modulus  and  pressure  derivative  of  the  bulk  modulus  are  summarized  in 
Table  3.  Figure  2  shows  a  plot  of  shock  compression  data  generated  by  these 
investigators  as  pressure  vs.  specific  volume.  In  this  figure,  the  data  points 
denoted  by  HEL  are  at  or  below  the  HEL;  PLl  indicate  inelastic  deformation  of 
SiC,  and  PL2  indicate  the  transformed  phase  of  SiC.  This  figure  shows  that  the 
compression  data  of  McQueen  et  al.  [10],  Marsh  [11],  and  Gust  et  al.  [12]  are 
consistent  wilh  each  other.  The  shock  compression  data  of  almost  fuUy  dense 
6H-SiC  reported  by  Sekine  and  Kobayashi  [13,  14]  is  stiffer  than  of  the  low- 
density  KT  SiC.  The  data  obtained  by  Sekine  and  Kobayashi  labeled  PL2  roughly 
follows  the  initial  trajectory  of  the  McQueen  et  al.  [10]  data  above  100  GPa.  The 
curve  in  this  figure  represents  the  compression  of  6H-SiC  obtained  by  using  the 
values  of  bulk  modulus  and  its  pressure  derivative  given  in  Table  3  and  relation 
(5). 

Table  3.  Density  (po),  HEL,  Bo,  and  B'ofor  SiCs  used  in  shock  compression  investigations. 


Type 

po 

(Mg/m3) 

HEL 

(GPa) 

Bo 

(GPa) 

B'o 

Pmax 

(GPa) 

Reference 

KT-SiC 

3.12 

8±3 

184-201 

2.8-3.4 

BKEEHi 

6H-SiC 

3.22 

18.0  ±0.7 

230^ 

4.6 

USSli 

3  Sekine  and  Kobayashi  [13, 14]  took  the  value  of  Bo  from  Bassett  et  al. 


[8]  for  their  analyses. 


Gust  et  al.  [12],  based  on  their  own  data  and  the  data  reported  by  McQueen  et  al. 
[10],  suggested  the  existence  of  two  phase  transitions  in  their  SiC.  The  first 
transition  in  SiC  was  estimated  to  occur  at  24  GPa  and  the  second  one  at  around 
96  GPa.  Sekine  and  Kobayashi  [13]  detected  only  one  phase  transition  in  6H-SiC 
at  105  ±  4  GPa  accompanied  with  a  volume  reduction  of  15  ±  3%.  The  transition 
was  not  complete  until  a  pressure  of  137  ±  4  GPa  was  reached.  The  structure  of 
the  new  phase  of  the  material  is  suggested  to  be  rock  salt  structure.  Further 
details  about  tiiis  transition  may  be  found  in  reference  [13]. 

2.3  Comparison  of  Hydrostatic  and  Hydrodynamic  Compression  of 
6H-SiC 


The  density  of  this  type  of  SiC  is  3.215  (Mg/ m^).  The  value  of  the  bulk  modulus 
(Bo)  at  the  ambient  condition  obtained  from  ultrasonic  wave  velocity 
measurements  is  221  GPa  (Table  2).  The  bulk  modulus  obtained  from  ultrasonic 
wave  velocity  measurements  is  the  adiabatic  modulus.  The  values  of  Bo  reported 
from  in-situ  x-ray  diffraction  measurements  at  static  high  pressure  are  the 
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Figure  2.  Shock  compressions  of  SiCs. 

isothermal  bulk  moduli  and  range  from  230  to  260  GPa  (Table  1).  The  ratio  of  the 
adiabatic  and  isothermal  bulk  modulus  of  SiC,  based  on  its  value  of  Gruneisen 
parameter  (1.2)  and  volumetric  thermal  expansion  coefficient  13.8  x  10-6  K-1  is 
1.005.  Thus,  the  difference  between  the  magnitude  of  bulk  modulus  obtained 
from  ultrasonic  wave  velocities  and  x-ray  diffraction  measurements  is 
insignificant.  No  measurement  of  bulk  modulus  is  available  from  the  shock 
wave  compression  data.  Sekine  and  Kobayashi  [13]  assumed  a  value  of  230  GPa 
for  the  bulk  modulus  reported  by  Bassett  et  al.  [8]  and  obtained  a  value  of  4.6  for 
B'o.  The  value  of  B'o  obtained  from  their  data  by  using  221  GPa  as  the  value  of 
bulk  modulus  determined  from  the  ultrasonic  velocity  measurements  is  4.8,  not 
very  different  from  the  value  of  4.6.  However,  a  regression  analysis  of  their  data 
indicates  that  the  uncertainty  in  the  value  of  B'o  is  ±0.4.  The  value  of  B'o  obtained 
from  the  hydrostatic  compression  measurements  is  2.9  ±  0.3  (Table  1).  It  is  of 
interest  to  note  that  the  estimated  values  of  B'o,  based  on  value  of  Gruneisen 
parameter  as  1.2  and  using  Slater  and  Dugdale-McDonalds  relations,  are  27  and 
3.4,  respectively.  Since  ultrasonic  wave  velocity  measurements  provide  the 
most  accurate  value  of  bulk  modulus  at  the  ambient  condition  and  are  more 
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likely  to  be  representative  of  6H-SiC  used  in  large-scale  applications  rather  than 
those  obtained  from  the  x-ray  diffraction  measurements  at  high  pressures,  it  is 
preferable  to  use  a  value  of  221  GPa  for  Bo  of  6H-SiC.  The  reason  for  this  is  that 
x-ray  diffraction  measurements  at  static  high  pressures  provide  the  compression 
of  SiC  itself  and  are  not  influenced  by  porosity  or  free  impurities  in  the  sample. 
A  larger  scatter  in  the  value  of  Bo  of  6H-SiC  reported  in  Table  1  may  arise  from  a 
host  of  sources,  including  the  one  related  to  the  differentiation  of  the  pressure 
volume  data.  Under  an  ideal  condition,  one  would  expect  the  value  of  the  bulk 
modulus  for  a  material  obtained  from  high-pressure  x-ray  diffraction 
measurements  to  be  consistent  with  the  ultrasonic  value  and  also  more  consistent 
with  one  another  than  reported  for  6H-SiC  by  Bassett  et  al.  [8]  and  Yoshida  et  al. 

[9]. 

The  compression  curves  of  6H-SiC  obtained  by  using  221  GPa  for  the  value  of  Bo 
and  2.9  ±  0.3  and  4.6  for  the  values  B'o  are  compared  with  the  compression  curve 
generated  from  the  results  of  compression  and  shear  wave  velocity 
measurements  to  characterize  the  shocked  state  of  6H-SiC  by  Yuan  et  al.  [20]. 
Figure  3  shows  the  results  such  calculations.  This  figure  shows  that  the 
compression  curves  of  6H-SiC  obtained  by  Yuan  et  al.  [20]  and  Sekine  and 
Kobayashi  [13]  are  almost  identical  to  one  another  to  60  GPa.  It  suggests  that  one 
may  use  the  equation  of  state  generated  from  the  data  of  either  of  these  two 
investigations  with  some  confidence  to  at  least  90  GPa.  Numerical  simulations 
are  generally  done  by  expressing  pressure  as  a  function  of  the  third  degree 
polynomial  of  p,  where  p  is  defined  as  {(Vo/V)-l)}.  Such  a  tnmcated  equation  of 
state  for  6H-SiC  using  the  data  from  Sekine  and  Kobayashi  [13]  may  be  given  by 

P  =  221  p  +  398  p2  +  283  p3.  (8) 

Figure  4  shows  that  relation  (8)  actually  represents  the  compression  of  6H-SiC 
very  well  to  60  GPa  and  satisfactorily  to  90  GPa.  This  suggests  that  relation 
(7)  can  be  used  to  calculate  compression  of  an  intact  6H-SiC  to  at  least  90  GPa. 

2.4  Shear  Strength  of  6H-SiC 

Shear  strength  of  SiC-B  under  shock  compression  through  simultaneous 
measurements  of  longitudinal  and  lateral  stress  measuring  has  been  measured 
by  Feng  et  al.  [19]  and  Bourne  et  al.  [18].  It  is  of  interest  to  calculate  the  values  of 
shear  stress  sustained  by  6FI-SiC  obtained  from  the  equation  of  state  given  by 
relation  (8)  and  the  longitudinal  stress  measured  near  and  above  its  HEL  in 
references  [16,  18,  19]  and  compare  them  witti  the  values  obtained  from  the 
measurements  of  lateral  and  longitudinal  stress  tmder  shock  loading.  The  values 
of  longitudinal  stress  and  the  associated  volume  change  needed  for  calculations 
were  given  only  in  the  paper  published  by  Feng  et  al.  [19].  The  relation  between 
the  longitudinal  stress  (ol),  pressure  (P),  and  shear  stress  (x)  at  a  given  strain  is 
given  by 

x  =  0.75{aL- P}.  (9) 
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Figure  4.  Comparison  of  compressions  of  6H-SiC  generated  from  relations  (5)  and  (8). 

Table  4  gives  the  results  of  the  calculations.  This  table  shows  that  the  calculated 
values  of  shear  stress  and  those  obtained  through  the  simultaneous 
measurements  of  longitudinal  and  lateral  stresses  agree  with  one  another.  The 
highest  magnitude  of  longitudinal  shock  experiments  were  done  by  Crawford 
and  reported  by  Feng  et  al.  [17].  Results  of  two  experiments  reported  in  their 
Table  3  are  given  in  Table  4  in  this  report.  The  values  associated  with  (a)  and  (b) 


10 


Table  4.  Measured  and  calculated  values  of  shear  stress  sustained  by  6H-SiC  under  plane 
shock  wave  compression. 


Experiment 

No. 

Stress 

_ (GPa) 

P- 

Calculated 

(GPa) 

Longitudinal 

Lateral 

Shear 

Pressure 

Shear 

Feng  et  al.  [19] 

1 

1.84 

4.18 

0.0209 

4.80 

4.05 

2 

12.90 

5.28 

0.0264 

6.12 

3 

5.80 

0.0312 

7.29 

4 

3.60 

6.20 

0.0336 

7.89 

6.09 

5 

5.10 

6.85 

0.0412 

9.80 

6.75 

6 

20.90 

6.94 

6.98 

0.0479 

11.53 

7.03 

7 

24.20 

10.40 

6.90 

0.0610 

15.03 

6.88 

Crawford's  experiments  reported  by  Feng  et  al.  [17] 

26.8 

0.0875 

22.8  ±  0.5 

3.0  ±0.4 

HiEISH 

39.6 

0.1314 

36.9  +  1.1 

2.0  ±0.8 

ma^ 

28.5 

0.0785 

20.2  ±  0.4 

6.2  ±0.3 

SC-4(b) 

40.8 

0.1237 

34.4  + 1.0 

4.8  ±0.8 

Grady  and  Kipp  [16] 

CE-4 

27.6 

0.075 

19.2  ±0.4 

6.3  ±  0.3 

CE-5 

36.3 

0.1087 

29.4  ±0.8 

5.2  ±0.6 

CE-Sl 

48.8 

0.1534 

Bourne  et  al.  [18] 

1 

16.7 

3.3 

6.7 

2 

21.2 

4.1 

8.6 

3 

23.4 

8.4 

for  these  experiments  are  based  on  the  calculation  of  longitudinal  stress  and  p 
using  die  elastic  perfectly  plastic  model  and  the  pressure-dependent  strength, 
stress  relaxation  model  proposed  by  Feng  et  al.  [17],  These  calculations  indicate 
that  die  values  of  shear  stresses  based  on  the  pressure-dependent  strength,  stress 
relaxation  model  are  consistent  with  the  other  values  of  the  shear  stress  given  in 
this  table. 

Since  the  experimental  results  of  Grady  and  Kipp  [16]  on  a-SiC  manufactured  by 
Eagle-Picher  was  shown  to  coincide  with  the  longitudinal  shock  of  SiC-B  in 
reference  [17],  the  values  of  shear  stress  sustained  by  this  material  are  also 
calculated  and  given  in  Table  4.  The  calculation  shows  a  gradual  reduction  in 
the  shear  strength  of  die  a-SiC  with  an  increase  in  the  longitudinal  stress.  Could 
this  be  related  to  the  observed  trend  in  the  change  in  ionicity  of  the  bonds  in  SiC 
and  the  development  of  covalency  in  SiC  at  around  40  GPa  by  liu  and  Vohra 
[31]?  The  measurements  by  Bourne  et  al.  [18]  appear  to  be  consistent  with  the 
measurements  of  Feng  et  al.  [19],  but  tend  to  be  a  bit  larger  than  their  values  at 


11 


longitudinal  stress  exceeding  21  GPa.  There  is  clearly  a  need  for  such 
measurements  at  and  above  20  GPa  in  SiC-B. 


3.  Concluding  Remarks  and  Future  Work 


The  findings  of  this  report  may  be  summarized  as  follows. 

•  The  equation  of  state  obtained  from  in-situ,  high-pressure,  x-ray  diffraction 
measurements  provides  very  precise  values  of  compression  (i.e.,  strain). 
The  somewhat  less  precise  measurement  of  pressure  may  be  due  to 
nonhydrostaticity  of  the  pressure  imder  which  such  measurements  are 
made.  A  very  good  example  are  the  measurements  of  compression 
obtained  by  Bassett  et  al.  [8]  and  Yoshida  et  al.  [9].  It  was  shown  that  while 
the  compression  loci  for  6H-SiC  obtained  in  these  investigations  were 
almost  identical,  the  values  of  bulk  modulus  and  its  pressure  derivative 
were  not  the  same  (Figure  1  and  Table  1).  Lack  of  knowledge  about  the 
complete  set  of  elastic  constants  for  a  single  crystal  of  6H-SiC  prevents 
making  any  judgment  about  the  quality  of  measurements. 

•  The  equation  of  state  obtained  from  the  shock  compression  measurements 
depends  on  the  porosity  and  free  impurities  in  materials  indicated  by  the 
measurements  on  KT-SiC  and  6H-SiC  (Table  3).  In  the  case  of  6H-SiC,  the 
investigations  by  Sekine  and  Kobayashi  [13]  and  Yuan  et  al.  [20]  showed 
that  the  compression  meastued  and  expressed  as  a  function  of  pressure  and 
mean  stress,  respectively,  are  identical,  and  the  value  of  bulk  modulus  used 
in  both  investigations  was  of  a  magnitude  similar  to  that  obtained  from 
ultrasonic  wave  velocity  measurements  (Table  2  and  Figure  3). 

•  The  compression  curves  obtained  from  in-situ,  high-pressure  x-ray 
diffraction  measurements  of  Bassett  et  al.  [8]  and  Sekine  and  Kobayashi  [13] 
are  consistent  with  one  another,  but  also  owe  their  consistency  to 
serendipity  in  the  sense  that  these  investigators  obtained  the  value  of  only 
one  of  the  two  parameters  (i.e..  Bo  and  B'o,  and  made  correct  assumptions 
for  the  value  of  the  remaining  one). 

•  The  experiments  of  Sekine  and  Kobayashi  [13]  yield  a  value  of  B'o  4.6  ±  0.4. 

•  The  third  degree  polynomial  equation  for  use  in  numerical  simulations  for 
6H-SiC  (i.e.,  SiC-B)  is  given  by 

P  =  221  p  +  398  p2  +  283  p3.  (10) 

•  The  procedure  used  to  estimate  shear  stress/strength  (x)  sustained  by  SiC-B 
from  hydrodynamic  compression  and  longitudinal  stress  measurements 
agree  with  the  values  obtained  from  the  measured  longitudinal  and  lateral 
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stress  measurements  (Table  4).  To  this  date,  the  procedure  has  always 
yielded  correct  results  for  other  materials  [32, 33].  The  value  of  t  appears  to 
be  constant,  even  though  the  collective  trend  indicates  that  SiC-B  may  be 
losing  its  shear  strength  with  an  increase  in  longitudinal  stress  beyond 
40  GPa. 

•  There  is  a  need  to  supplement  the  existing  measurements  of  shear  strength 
of  SiC-B  to  vahdate  the  gradual  reduction  of  shear  strength  at  shock 
stresses  beyond  24  GPa. 

•  There  is  no  experimental  data  on  any  type  of  SiC  to  provide  information 
pertaining  to  the  effect  of  shock-induced  damage  on  its  subsequent 
dynamic  response.  Shock-reshock  experiments  on  SiC,  specifically  SiC-B, 
wiU  complement  the  existing  data  [34, 35]. 

•  The  extensive  information  available  about  the  spall  behavior  of  SiC-N 
needs  to  be  supplemented  by  investigating  shock  compression  and  release 
response  of  SiC-N  above  its  HEL. 

•  Finally,  nonconvergence  of  tire  shock  compression  data  for  KT-SiC  and 
6H-SiC  in  Figure  2  suggests  that  it  may  be  worth  examining  and  analyzing 
ihe  shock  compression  data  of  McQueen  et  al.  [10]  and  Marsh  [11]  on  a-SiC 
of  densities  2.9, 3.0,  and  3.1  Mg/m^  to  look  into  the  feasibility  of  using  their 
data  to  obtain  compression  response  of  damaged  SiC.  The  assumption  is 
that  the  different  densities  of  SiC  represent  differing  amounts  of  damage  in 
the  material. 
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